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We show that a vortex current is created around a skyrmion spin texture under magnetic field 
due to a radial spin motive force in a two-dimensional metal with localized magnetic moments even 
in the absence of any superconductivity correlations. The effect is expected both for ferromagnetic 
and for antiferromagnetic systems. The Skyrmion-induced vortex mechanism provides a picture for 
large Nernst signals observed in the pseudogap phase of the high- Tc cuprates. 
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PACS numbers: 72.25.-b,72.15.Gd,74.72.-h 



Skyrmions were originally proposed to describe 
baryons in terms of meson fields [l|. The central idea 
of constructing a topological object starting from an un- 
derlying field theory is widely used in many areas of 
physics [2[ . In condensed matter physics the existence of 
skyrmions is suggested and confirmed inpart in quantum 
Hall systems with filling fraction ^ ~ 1 [3|,|j], spinor Bose- 
Einstein condensate JS, itinerant ferromagnets [6J, ne- 
matic liquid crystals [Tj , and high-temperature supercon- 
ductors [si |9| . From the topological nature of skyrmions 
the Berry phase effect [10[ leads to an effective mag- 
netic field, which has been studied intensively in vari- 
ous systems llj . In this Letter I show that under mag- 
netic field skyrmions create a vortex current due to spin 
motive force without relying on any superconductivity 
correlations. The skyrmion-induced vortex mechanism is 
applicable not only to ferromagnetic systems but also to 
antiferromagnetic systems. 

In a two-dimensional isotropic ferromagnet, where the 
local magnetization direction vector is represented by 
m = (toi, TO2, 'tis), a skyrmion spin texture at the origin 
has the following form. 
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where r — \J x^ -|- y^ and A is the core size of the 
skyrmion. At infinity m approaches to the magnetization 
direction, which is taken as 63. (Hereafter Cj represents 
the unit vector in the j'-axis.) This skyrmion spin texture 
is characterized by the topological charge, 

Q-^ /"d^rm(r)- [a^m(r) xaj,m(r)]. (2) 

If conduction electrons interact with m through an ex- 
change coupling 



He - MJ, f Sv [V'^ (r) CT^ (r)] • m (r) , 



(3) 



where M is the magnetization and the vector cr is, cr — 
(cTi , cr2 , tTa ) with Oj the Pauli matrices, the Berry phase 
effect associated with the non-vanishing scalar chiral- 
ity density m (r) • [Qa^m (r) x dyUv (r)] produces an ef- 
fective magnetic field effect for the conduction electrons. 



Anomalous Hall effect arising from such a Berry phase 
effect has been discussed in manganites [l2|, pyrochlores 
[13], and MnSi[i3. 

In this Letter, we study another Berry phase effect 
that emerges under magnetic field. The effect is associ- 
ated with the spin motive force discussed in the context 
of spintronics 15[ and experimentally observed recently 
[l6| . If there is a domain wall, the spin motive force is in- 
duced along the domain wall under magnetic field. Here 
the same analysis is applied to skyrmion spin textures. 

We start with a ferromagnetic metal. 






v)K{-inV)i^s{r)+'H, + 'H^, (4) 



where K is the kinetic energy operator. The term T-Lra de- 
scribes the interaction for m, and contains Dzyloshinsky- 
Moriya interactions. Here we assume that ■Hm stabilizes 
a skyrmion texture for m, and do not discuss explicit 
mechanisms for stabilizing a skrymion 17 1. 

The coupling term T-Lc tends to align the conduction 
electron spin to m. In order to include this correla- 
tion effect, we perform an SU(2) gauge transformation, 
i}{T) -^ C/(r)V'(r), where [/ (r) = n (r) • cr [13, [I^] with 



n(r) 



(mi (r) , m2 (r) , 1 + TO3 (r)) 
y/2 [1 + mg (r)] 



(5) 



The gauge potential associated with the SU(2) gauge 
transformation is given by a = —iTiU^'VU and uq — 
ihU^dtU . In terms of n, we obtain 



a = ?i(T • (n X Vn) . 



ao 



-ha ■ (n X dtn) 



(6) 



(7) 



Note that 9tn is computed from the Heisenberg equation 
of motion. From these expressions we find the SU(2) 
gauge field. 



-2h(T ■ {dtn X Vn) . 



b = fiEijkai {Vuj) X (Vnfe) , 



(8) 



(9) 



where Sijk is the antisymmetric tensor. 

Now we calculate the SU(2) gauge field created by a 
skrymion at the origin, Eq. ([T}. Introducing the cylin- 
drical coordinate (r, 0, z), we obtain 
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and e = 0. It is easy to see that bzu = [&z]tt is 
equal to m • {dxm x dyva) /2 [12j. The field &z„ be- 
haves like a magnetic field for the conduction elec- 
trons. For A = lOOA, the average field strength is 
J^dr2Trrbzu/{2TTe\'^) ~ 3.3 T. With decreasing A, b,u 
increases, and so bzu can be very large. However, this 
effective field is not directly observable because it is not 
a real magnetic field. 

Now we apply the magnetic field B = (0, 0, —B) to the 
system. From the Zeeman energy term the dynamics of 
m is given by 



d_ 
dt 



m — aB X m, 



(11) 



where a — gMfiB/fi with g the g-factor and fiB the Bohr 
magneton. Using the solution of this equation of motion, 
we find 



n(r,i) 



A cos $ (i) ei + A sin $ (i) 62 



res 



\fr^ 



A2 



(12) 



with $ (i) = (j)—aBt + (j)o. {(f)Q is a constant.) The field b 
is given by Eq. PT!]) with being replaced by $(i). Now 
the field e is nonzero and given by 



e = 



2aBh\^ 

(r2 + X^f 
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(13) 



Note that the effective field e is created in the radial 
direction [20J. This is understood as follows. If one sees 
the skyrmion spin texture along a straight line crossing 
the origin in the x-y plane, the spin configuration looks 
like a domain wall. From the analysis of the domain 
wall we know that the spin motive force is created along 
The same is true in each direction. Thus, the 
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spin motive force e is created in the radial direction for 

skyrmion spin textures. 

Since e^u is non-zero and there is the external magnetic 
field and bzu, a drift of conduction electrons is induced. 
The drift is circular around the skyrmion because ezu is 
in the radial direction. Thus, the drift motion leads to 
the vortex current. The drift velocity is 



v^ 



2aB\^i 



2A2 + {eB/h) (r2 + A^)" 



(14) 



Figure [T] shows this drift velocity as a function of r/A. 
For r/A < 1, bzu plays a major role for the drift. While 
for r/A > 1, B plays a major role. For B > Q the ex- 
ternal magnetic field is parallel to bzu- If we consider 
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FIG. 1. (color online) The drift velocity as a function of r/A 
for different magnetic fields at A = 40^4 with vq — gMc with 
c the speed of light. The inset shows the drift velocity versus 
r/A for different A (in units of A ) at S = 10 T. 



an auti-skyrmion, the drift velocity changes its direction 

at r — J \ [V2iB — a) with Ib — \/eB/h the magnetic 
length. We do not consider this case in the following anal- 
ysis because such a current distribution is energetically 
unfavorable, though it is not forbidden. 

The circular drift velocity (1141) leads to a vortex cur- 
rent. The magnetic field By created by this vortex cur- 
rent is computed as in that for vortex current in a layered 
superconductor 21]: 



B„ 



dqqjQ{qp)Ao{q)ex.p{-q\z\), (15) 



^o(g) = -- 



27r 



dppJi [qp) (J2d)^ , (16) 



with J„(x) the Bessel function and the two-dimensional 
current density is (J2d)j, = e.n2dV^. Here n2d is the two- 
dimensional conduction electron density. This magnetic 
field By is shown in Fig. [2j The field is measured in units 
of tcsla and is normalized by a dimensionless parame- 
ter Bo = gMn2dA2. Note that usually vortex currents 
are induced if there is a superconductivity correlation. 
However, the skyrmion-induced vortex current does not 
require any superconductivity correlations. 

The same formula can be applied to the antiferromag- 
netic case as well. Let us consider a square lattice and 
introduce A and B sublattice. For this case we take m(r) 
as the local staggered magnetization direction vector. 
The conduction electron spin residing at Rj is aligned 
to m(Rj) by the unitary transformation matrix [/(Rj). 
Carrying out this transformation, we obtain the same 
Hamiltonian except for the exchange interaction term. 
The sign of the exchange interaction term is different for 
A and B sublattice. In the strong coupling limit, the 
conduction electron spin at B sublattice is anti-parallel 
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FIG. 2. (color online) The magnetic field Bv created by 
the vortex current for difTerent magnetic field at z = with 
A = 40A. The inset shows B„ for A = 10, 20, 40, 60 (in units 
of A) from top to bottom. 



to that at A sublattice. The relevant components of the 
SU(2) fields are 



[e]ti = ([e] It)* = 



2aBhX^ 

[r^ + X^f 
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(17) 



(18) 



Note that the SU(2) gauge fields acting on the conduction 
electrons are not staggered, in spite of the fact that the 
magnetization direction vector is staggered. Denoting 
e = ^ ■ Bjaj and b = ^ bjCj, we find the gauge field 
components e^ and b^. For a bond vector connecting 
from B sublattice to A sublattice, these components take 
opposite sign compared to the case of a bond vector con- 
necting from A sublattice to B sublattice. A crucial point 
is that this sign change cancel completely in the drift cre- 
ated by e and b. The situation is clearer if we carry out 
an additional Unitary transformation V = exp (i^crz/2) 
that makes e 



ti 



e^t and b^; 



h^-f and removes the 



sign change. The drift velocity calculation is similar to 
the ferromagnetic case, and we obtain 



2aBX^ 



Vff, 



2X^ + {eB/n)r{r^ + X^y 



(19) 



The difference arises from the relevant component of the 
gauge fields as stated above. Figure [3] shows the mag- 
netic field created by the drift calculated by Eqs. (fT5|) 
and pB|) using Eq. ([TO| . Although r dependence of the ef- 
fective fields are different between the antiferromagnetic 
case and the ferromagnetic case, the results are similar. 
For r < X, the effective field b plays a major role for 
the drift. In this regime, the drift velocity has the same 
form as that for the ferromagnetic case. This is the rea- 
son why we obtained similar results for both cases for 



FIG. 3. (color online) The magnetic field _B„ as a function 
of r/X for the antiferromagnetic case. The value of B^ is 
computed for different magnetic field at z — with A = 40 A. 
The inset shows i3„ for A = 10, 20, 40, 60 (in units of A) from 
top to bottom. 



r < X. Meanwhile for r > X the effective electric field is 
suppressed by the factor X/r. Thus, By decreases rapidly 
compared to the ferromagnetic case. 

In order to stabilize the vortex current, we need to 
keep the radial effective electric field e. However, in 
general a radial electric field is screened by conduction 
electrons. Therefore, to stabilize the vortex current the 
Thomas-Fermi length, which is proportional to l/n2dj 
should be larger than A, otherwise the effective electric 
field is screened by the conduction electrons. 

Now it is natural to ask a question: What is the 
difference between the vortex current arising from the 
skyrmion spin texture and a conventional vortex aris- 
ing from a superconductivity correlation? As long as 
the external magnetic field is kept, the skyrmion- induced 
vortex current continues to flow even in the presence of 
dissipations. However, the current decays rapidly when 
one turns off the magnetic field. This makes the crucial 
difference compared to conventional vortices in a super- 
conductor. In addition, the flux quantization occurs only 
for a mesoscopic sample. 

In the formulation above we have assumed that there 
is magnetic long-range order. But this assumption is not 
necessary. Suppose the system does not have magnetic 
long-range order and let ^ be the magnetic correlation 
length. For ^ 3> A, we may apply the formula above to 
the system as well. If £, is comparable to A, the field 
b and e are multiplied by the factor exp(— 2r/^), and 
M should be replaced by the average of spins over the 
domain. The amplitude of the fields would decrease by 
these factors but one may still expect a vortex formation. 

So far we have assumed that skyrmions are static. 
However, it is possible to consider a moving skyrmion 
[22] and compute a vortex current around it. For a mov- 
ing skyrmion, the antiferromagnetic case is much better 
than the ferromagnetic case. If one adopts the non-linear 



a model for the description of the skyrmion spin texture, 
the skyrmion is a soliton solution of the non-linear u 
model. However, for the ferromagnetic case the skyrmion 
spin texture is stable only for the static case. Because of 
the quadratic dispersion of spin waves, the chiral nature 
of the skyrmion spin texture is lost upon skyrmion prop- 
agation. By contrast, for the antiferromagnetic case a 
moving skyrmion spin texture solution is constructed by 
a Lorentz boost. This is possible because the non- linear 
(J model has a relativistic form. Of course there are devi- 
ations from the relativistic dynamics in the real system, 
and those deviations would lead to a finite life time even 
in the antiferromagnetic case. 

In order to observe the vortex currents predicted above 
candidate systems are Fei_a;Co2;Si, where real-space ob- 
servation of a two-dimensional skyrmion crystal was re- 
ported [23|, and MnSi, where magnetoresistance mea- 
surements suggest the presence of skyrmions |14| . An- 



other interesting application of the theory is that for 
cuprate high-temperature superconductors. In recent ex- 
periments large Nernst effects were observed [24] above 
the superconductivity transition temperature T^- For 
the single band system, which is believed to be the case 
for the high-Tc, the Nernst effect is absent because of 
Sondheimer cancellation [2J|. One scenario for the large 
Nernst effect above T^ is to assume the presence of pre- 
formed Cooper pairs. Assuming a skyrmion texture car- 
ried by a doped hole provides an alternative scenario. To 
make clear whether the skyrmion spin texture is formed 
in the high- Tc cuprates or not, the best way would be 
to investigate Li-doped underdoped samples. Establish- 
ing a skyrmion spin texture by observing the vortex cur- 
rent around a doped hole bound to a Li ion provides a 
crucial step to uncover the mechanism of cuprate high- 
temperature superconductivity [25|. 

To conclude, we have discussed a mechanism of vor- 
tex current formation under magnetic field based on the 
skyrmion spin texture without relying on any supercon- 
ductivity correlations. The vortex current is created for 
antiferromagnetic systems as well as ferromagnetic sys- 
tems. The presence of the vortex can be verified by ob- 
serving the magnetic field created by the vortex current. 
Searching for a skyrmion-induced vortex current in un- 
derdoped cuprates would be a key step to justify the 
relevance of the skyrmion spin texture in high- Tc. 
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